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Abstract: Following a suggestion by Smit, the CP odd terms of the effective action of 
the Standard Model, obtained by integration of quarks and leptons, are computed to sixth 
order within a strict covariant derivative expansion approach. No other approximations 
are made. The final result so derived includes all Standard Model gauge fields and Higgs. 
Remarkably, at the order considered in this work, all parity violating contributions turn 
out to be zero. Non vanishing CP violating terms are obtained in the C-odd P-even 
sector. These are several orders of magnitude larger than perturbative estimates. Various 
unexpected regularities in the final result are noted. 
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1. Introduction 



Whereas CPT symmetry is preserved of necessity in any theory with Lorentz invariance, 
locaUty and unitarity Q no such mechanism is at work to preserve CP Q. Nevertheless, 
CP is very weakly broken in the Standard Model ^. No breaking is detected in the 
strong sector, where the coupling constant of the operator G^yG^" is compatible with zero 
Q. In the electroweak sector no breaking is possible for less than three generations and 
even in this case the breaking would not occur in the presence of mass degeneration, as 
is the case in the lepton sector for massless neutrinos So there is CP violation in the 
Standard Model but it is rather small, as compared to the maximal breaking of C or P in 
the electroweak sector. 

The amount of CP breaking is relevant in early universe baryogenesis ||^]. Baryon 
asymmetry generation is assumed to take place near the electroweak transition tempera- 
ture, T of the order of 100 GeV. At such high temperatures quark masses, except that of 
the top quark, can be considered small. This suggests to treat them as a perturbation. 
Since these masses follow from the Yukawa coupling of quarks to the Higgs field, this is 
equivalent to treat those Yukawa vertices perturbatively. However, as noted above, CP 
violation is elusive as no CP breaking term can be produced at low orders. The simplest 
such term is the Jarlskog determinant which appears at order twelve 

A = J(m^ - ml){ml - m^){mj - ml){mj - m^)(m^ - mf)(mf - mj) . (1.1) 

where J is the dimensionless Jarlskog invariant, constructed with the Cabibbo-Kobayashi 
Maskawa matrix. This term has dimension twelve and by dimensional counting it is usually 
assumed to enter in the effective action scaled by T^^. At the electroweak transition 
temperature the ratio A/T^^ is extremely small, of the order of 10~^^. If this estimate is 
correct this poses a problem to account for baryogenesis using the Standard Model |^. 

A simpleminded transcription of the above estimate to the zero temperature case can 
be achieved by simply replacing the scale T by the Higgs field condensate v = 246 GeV, 
A/v^"^ = 10~^^. This is several orders of magnitude smaller than CP violation as measured 



in meson decays, where dimensionless parameters are of the order of 10~^ |10|. However, 
the direct use of the Jarlskog determinant is not justified at zero temperature, where the 
quark masses can no longer be treated perturbatively. 



Smit |11| made the observation that a non a perturbative treatment would yield in a 
natural way much larger couplings for CP breaking operators, as such couplings would come 
out as rational functions (with logarithms) of the quark masses times the Jarlskog invariant. 
Specifically it was proposed to study the effective action of the Standard Model obtained 
after integration of the fermions in the theory. The full functional is, of course, beyond 
an exact computation and some type of classification and selection of the resulting terms 
is required. The proposal was then to organize the terms within a covariant derivative 
expansion, which being non perturbative has the potential of yielding a more reliable 
estimate for the couplings. 

The study of the leading order terms in the abnormal parity sector was undertaken 
in 1 11]. This is the sector driven by the Levi-Civita pseudo-tensor and includes the Wess- 
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Zumino-Witten term. In the Standard Model this is equivalent to the parity odd sector. 
The corresponding operators have dimension 4, counting only the dimension carried by the 
dynamical fields except the Higgs. Unfortunately no non vanishing contribution was found 
to fourth order, although CP breaking contributions were expected at dimension 6 |11|. 
Such dimension 6 operators, with non vanishing coupling, have been found in [12| where 
also the abnormal parity sector was studied. As expected, the CP violating term found is 
indeed sizable as compared to perturbative estimates. 

The two calculations of CP violating terms in the Standard Model effective action 
just described are based on the technique introduced in [13, 14]. In this approach the full 
effective action for a generic theory of fermions coupled to chiral gauge fields is computed 
within a strict derivative expansion, to a given order. |11| uses directly the result in |14| 
which holds to fourth order, and particularizes it to the Standard Model while selecting 
just the CP breaking terms. |12| carries out the same reduction from general to particular 
after extending the generic calculation to sixth order. This is done using the worldline 
formalism to deal with Dirac traces and momentum integration, as explained in |15|, as an 
alternative to do the same thing with the techniques applied in |14|. 

In the present work we also undertake the calculation of the CP odd component of the 
Standard Model effective action at zero temperature derived by integration of quarks and 
leptons, and organize the terms so obtained by means of a covariant derivative expansion. 
However, unlike previous calculations, ours is carried out from scratch by applying the 
recently derived technique described in |]l^, |l^]. The difference with previous approaches 
is that we particularize very early our treatment to the Standard Model and this allows us 
to select from the beginning terms which are candidate to break CP invariance and neglect 
irrelevant CP even terms. This is useful as CP breaking imposes very restrictive conditions 
and selects very few candidates. Another difference is that we consider terms of normal 
parity (i.e., P even and consequently C odd) as well as of abnormal parity (P odd, C even). 
Seemingly, in the literature, the CP breaking terms have been assumed to be of abnormal 
parity only. Perhaps this is because in the so called strong CP problem, the CP breaking 
terms involve the Levi-Civita pseudo-tensor in the gluon sector or 75 in the quark sector 
lH, H]. Mor eover, the simplest CP odd terms one can write for the effective action are 
also odd under intrinsic parity, e.g., ii {F^yF^^^) ||2^, where F^^, is constructed with the 
covariant derivative of the electroweak group SU(2) x U(l). In addition, the CP odd and P 
odd sector is of interest in the study of electric dipole moment. Nevertheless, it is perfectly 
possible to write down operators, constructed with the Standard Model gauge fields and 
the Higgs, having normal intrinsic parity and odd under CP transformations. 

We carry out a detailed calculation including all Standard Model fields and all CP vi- 
olating contributions to order six in the covariant derivative expansion. Diagrammatically, 
this consists of all one-loop Feynman graphs with fermions running on the loop and up to 
six gauge fields or derivatives (four-momenta) attached as external legs and any number of 
Higgs fields coupled to the quark or leptons. The result is given in the unitary gauge. The 
calculation presented is largely self-contained and some of the main conclusions obtained 
can be reached by a by-hand computation. This is the case for our most unexpected find- 
ing, namely, that there are CP violating terms of order six in the normal parity sector of 
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the Standard Model, but all terms vanish in the abnormal parity sector. This is at variance 
with the result in [l^. It is not clear to us from where the discrepancy arises as the two 
calculations are conceptually similar although technically different. In any case we have 
double-checked our results to confirm this conclusion. 

The paper is organized as follows. Section |^ summarizes a number of definitions, for- 
mulas and techniques relative to generic chiral gauge theories which will then be applied 
to the Standard Model. In section ^ we cast the fermionic sector of the Standard Model in 
the format previously described for generic chiral gauge theories. In section |^ features of 
CP violation at the level of the effective action are discussed. In section ^ the relevant mo- 
mentum integrals which appear in the calculation, as well as the associated selection rules, 
are obtained. Section ^ discusses the covariant derivative expansion within our approach 
to the effective action of the Standard Model. In section ^ the chiral invariant approach 
devised for generic theories is applied to the Standard Model in a way that allows to easily 
remove irrelevant CP even terms. Section ^ reproduces the previous result in |11] verifying 
that there are no CP breaking terms driven by operators of order four in the abnormal 
parity sector, and this result is extended to the normal parity sector as well. In section ^ we 
present a preliminary calculation for the particular case where Higgs field derivatives are 
neglected. This is of interest as this covers the result obtained in |12]. The cancellation of 
the abnormal parity contribution is made manifest there within a transparent calculation. 
Section IC presents the full result of our computation, Eq. ( 10.2 ). Various surprising regu- 
larities in the result are noted. The form of the loop function controlling the CP breaking 
operator is discussed in section 11 and how it is affected by infrared enhancement in the 



physically relevant chiral limit. Our conclusions are summarized in section 12 



2. Chiral gauge fermions 

In this section we collect important practical results relative to generic chiral gauge fermion 
theories which will be applied subsequently to the Standard Model. 

2.1 The effective action 

The Lagrangian describing the coupling of spin 1/2 fermions ('0) to chiral gauge fields 
{^l,r) and spin zero fields {rriLR, mRi) can be cast in the general form 

C{x) = i){x)i'D{x)tp{x), (2.1) 

where^ 



iD(x) = {i^ -y^{x))PR + Vl{x))Pl - mLR{x)PR - mRL{x)PL, (2.2) 

and 

Pr = \{1+1^), Pi = 1(1 -75), (2.3) 



"""We use Minkowskian signature (H ) and {7'', 7"} = 2(7^", 75 — , £0123 = +1. 
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and the external fields itilr, itirl, Vr,l are matrices in some internal space. Unitarity 
requires 

m{j,{x)=mRL{x), vIJx) = Vr,^{x), vIJx)=VlA^). (2.4) 
It will prove convenient to write the Lagrangian in matricial form, namely, 

£W = {fc,fe) f. . (2.5) 

Y^-Vr{x) -rriRL J \iPlJ 

This form is useful to expose the action of chiral transformations. 
The integration of the fermions provides the effective action 

Z = e^^ = J D^Dil) ^ '^"'^ ^^^^ ^^^^ , (2.6) 

which (modulo UV ambiguities) is given by 

ir[mLR,mRL,VR,VL] = TrlogiD. (2.7) 

The effective action can be decomposed into a normal parity component, T'^ (without 
Levi-Civita pseudo-tensor) and an abnormal parity component T~ (with the Levi-Civita 
pseudo-tensor) : 

r = r+-Fr-. (2.8) 

are also even/odd, respectively, under the pseudo-parity transformation, which can be 
defined as the exchange of the labels LR in the external fields. 

Among other symmetries, the effective action is invariant under the transformation 

rriLRix) mlji{x), mRL{x) m%L{x), 

VrA^) ^ -7r/V^,,(x), VlAo^) ^ -7r/li%(x), (2.9) 

where vr'^jy = diag(l, — 1, — 1, — 1) and x = {x^,—x) = ttx. This represents a CP trans- 
formation which we shall call full CP transformation to distinguish it from the physical 
one (which only acts on dynamical fields and not on parameters of the Lagrangian) ? The 
important property of this transformation is that it does not mix different chiral sectors. 

2.2 Euclidean space 

For convenience we shall work in Euclidean space, reverting to Minkowskian space at the 
end. 

^Likewise, full parity is also a symmetry of the effective action. It consists of exchanging LR labels and 

simultaneously x ^ x. Pseudo-parity, which can be defined as any of these two transformations without 
the other, is not a symmetry for general background field configurations. 
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The Euclidean metric is 6^^, so we put Euclidean indices as subindices. The pass from 
Minkowskian to Euclidean variables is achieved by the replacements^ 

■ip{x) ip{x), 'ip{x) — > ipix), 
niLRix) rriLBix), mRL{x) niRL^x), 
VRfiix) Vr^o{x), Vl^o{x) Vl^o{x), 
VR,i{x) -iVR^iix), VL,iix) -iVL,i{x), 

7° ^70, 7*^^7i; 75^ 75 = 70717273- (2.10) 
These replacements imply C{x) — > —^{x) and so 

= J D^jDiJje'-f'^"^^ ^ e"^ = g-/'^''^^^^) (2.11) 



with 

/:(x) = ^(x)DV(x) , (2.12) 

and 

D = + fi? ix))PR + {^ + yL {x))Pl + mLR{x)PR + mRL{x)PL. (2.13) 
Also, in Euclidean space 

V[mLB.mRL,VR,VL]=TY\og'D. (2.14) 
With the above prescriptions, in Euclidean space unitarity becomes 

m^j^ix) = mRLix), V^x) = -VrJx), V^x) = -VlA^), (2.15) 
while the full CP transformation becomes 

rriLRix) mlR{x), mRiix) mjji(x), Vr^x) vr^i.VS,^(x), VL,^ix) tt^^VI^x), 

(2.16) 

with TT^i, = diag(l, —1, —1, —1) and x = {x^, —x) = nx. 

An important property of the effective action in Euclidean space is that the normal 
parity component, r+, is real, and the abnormal parity component, T~ , is purely imaginary 
1 21]. This property is a consequence of unitarity and holds at the non perturbative level. 

2.3 Chiral invariant approach to the effective action 

The (Euclidean) Lagrangian ( 2.12| ) is invariant under local chiral transformations: 

Si^V V^'h ^bl) yo "i/ ^ 

However, as is well known, the corresponding effective action T[mLR, turl, Vr, Vl] displays 
an anomalous variation under chiral transformations |^, 23, 24 1. The anomaly has a 



^This corresponds to r;4 = 1 in Ref. ]l4 
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universal, geometrical, form and is saturated by the gauged Wess-Zumino-Witten (WZW) 



term, rgwzw["iLiJ) "^RL; ^l]; which also has a known geometrical form |25, This 
means that if the WZW term is subtracted from the effective action, the remainder, Fc, is 
chiral invariant: 

r = Fgwzw + Tc- (2-18) 
Let us note that Fgwzw contributes only to the abnormal parity component. Therefore, 

r+ = r+. 

Recently, it has been shown that the remainder Tc can be expressed in a convenient 
form in which chiral invar iance is manifest [^]. Indeed, let 

K = miRmRL- mj^l If)^ turl, (2.19) 

where -Dl,/^ = 9^ + and -Di?,/^ = + Vr^^, and so if is a second order differential op- 
erator. (Note that itT'rl ''^RL stands for the product of four consecutive operators.) 
Then"^ 

r+ = -^ReTr(logK), 

r- = -^iImTr(75logK) + rgwzw. (2.20) 
2.4 The derivative expansion 



The functional trace in ( 2.14 ) can not be computed in closed form in general. This suggests 
to use instead some systematic expansion to address its determination. In the derivative 
expansion scheme the effective action contributions are classified according to the number 
of covariant derivatives they carry. In this counting the spin zero fields, rriLR, turl, count 
as order zero whereas the derivatives or gauge fields, Vl,^, Vr^^, count as first order. 
Technically, this means to consider the family of Dirac operators 

T>t = (^ +t YR{tx))PR + (^ +t Yl {tx))PL + mLR{tx)PR + mRL{tx)PL, (2.21) 

and then expand the corresponding effective action in powers of the parameter t. After 
extraction of a global factor l/f^ in d space-time dimensions, the terms of order n contain 
n covariant derivatives (or gauge fields). At the diagrammatic level, a term of order n 
represents a one-loop Feynman graph with any number of scalar fields and n gauge fields 
as external legs, all of them at zero four-momentum, or less gauge fields as external legs 
and correspondingly more powers of the external momenta. We emphasize that gauge 
fields are assimilated to derivatives in such a way that they are of the same order. This 
ensures preservation of gauge invariance and as a consequence each order of the derivative 
expansion of the effective action is separately invariant under gauge transformations. 

Several remarks can be made: i) In even-dimensional space-times and at zero temper- 
ature, the derivative expansion of the effective action contains even orders only, ii) The 
abnormal parity sector starts at order d m. d dimensions, since it contains the Levi-Civita 



*In the notation of jl^], K here is Kl there, and we have used that Kn and K]^ are related by a similarity 
transformation. 
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pseudo-tensor, iii) In d dimensions, the terms beyond order d are ultraviolet (UV) con- 
vergent, and so free from UV ambiguities introduced by the renormalization. iv) Tgwzw 
contains only terms of order d (in d dimensions). I.e., within the derivative expansion, the 
WZW term vanishes at all orders beyond the lowest order one. This means that = T~ 
except at lowest order. And v) within the derivative expansion, the Minkowskian version 
of the effective action is real in the normal parity and in the abnormal parity sectors. The 
derivative expansion is an expansion around small external four-momentum, and so it can 
not reach the analytical cuts related to particle production. 

2.5 The method of symbols 

A convenient technique to carry out calculations within the derivative expansion is the 



method of symbols |27, 28 1. The method has been extended to curved space-time |29, 30 1 



and finite temperature [31, 32 1. 



This method will be used below in our calculation of the effective action. 

For any pseudo-differential operator / of the form f{D, M), constructed with covariant 
derivatives and external fields M{x) (all this non abelian in general) the method of 
symbols states, for the diagonal matrix elements of /, 



= / ^{x\f{D + p,M)\0) (2.22) 



{d is the dimension of the space of x). |0) represents the constant state (x|0) = 1, so that 
d^\0) = 0. For notational convenience, here and in what follows, we use a purely imaginary 
momentum = ik^ {k^ real), but denotes —PfiPfi = k"^ and d'^p = d'^k. The matrix 
element {x\f{D M)|0) just coincides with the standard symbol of /, as defined in the 



theory of pseudo-differential operators |33]. 



Invariance under the shift p^ ^ Pf_i + in the momentum integral ( 2.22| ) implies that 



can appear only in the form [D^, ] in the final integrated expression. This ensures gauge 
covariance of the right-hand side of ( 2.22| ), consistently with the obvious gauge covariance 



of the left-hand side. 

From ( ^.22 ) the derivative expansion is easily obtained just by formally expanding in 



powers of D^.^ The final step is to move all derivatives to the right (derivating everything in 
passing as dictated by Leibniz's rule) and to verify that terms with derivatives at the right 
(which would break gauge invariance) vanish after carrying out the momentum integration. 
The method of symbols then provides the (functional) trace of / as 

^(/) = / ^j0?^'f(D + P,M). (2.23) 
(The brackets {x\ |0) are usually omitted.) 



^p^ is the momentum running in the quantum loop and is of order zero in the derivative expansion. On 
the other hand Dp now acts only on the other external fields present in / and counts as first order. Direct 
expansion in powers of in f{D, M) (i.e. before applying the method of symbols) would not produce UV 
convergent integrals and would not correspond to the derivative expansion, as -D,i would contain not only 
the momenta of the external fields but also that of the quantum field running in the loop. 
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The method of symbols just described is simple enough but has the drawback that 
gauge invariance is not manifest during the calculation. An alternative technique, which 
we shall also employ in this work, is that of covariant symbols |3^, which has the 
virtue of being manifestly covariant from the beginning. Quite simply, the method consists 
in applying a similarity transformation in ( p.22 ), which changes nothing: 



{x\f\x) 



I 



{2'kY 

This can be seen to be equivalent to 

{x\f\x) 



f{D + p,M)e 



(27r)a 



f{D,M), 



where Du and M are manifestly gauge covariant: 



M 



(2.24) 



(2.25) 



2!^""'"'^^ap« 3l'^''''^''^^"dp^dpf, 
and these series are truncated at the desired order in the derivative expansion. 



(2.26) 



3. Standard Model 



3.1 Fermion sector of the Standard Model 

We can apply the previous general results to the Standard Model (SM) for quarks and 
leptons coupled to gauge fields and Higgs p^, Il0|. For quarks (and in Minkowski space) 



^SM,q(2;) = q{x) iDsM.q qix) . 



(3.1) 



In the notation of (|2.5|) the quark field, q{x), is a column matrix in the space Ti.LR ' 



n 



ud 



H 



gen 



has dimension four. 



"^coior ^Diraci where TClr = Ti-L ® Hr (chirality) has dimension two, 
'Hud = W„ © TLd (u or d quark type) has dimension two, Tigen = Ti-i © 7i2 © Ti.^ (generation) 
has dimension three, Wcoior has dimension three and TCbi 

( 

q{x) = 



q{x) = [uL,dL,UR,dR 



(3.2) 



dR 

UL 

Here only the LR and ud spaces are explicit while generation, color and Dirac indices have 
been left implicit. Likewise, in the unitary gauge, the Dirac operator takes the form 



^DsM,q = 





-^,4>Yd 



(?-\g%-y$-bs^a$a 






l_ 

(3.3) 
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Yuxi are 3x3 matrices in generation space which denote the Yukawa couphngs of the 
quarks with the Higgs field ^(x). Ga,^ are the gluon fields, with coupling constant gs, and 
Aa the Gell-Mann matrices in color space. is the U(l) weak hypercharge gauge field, 
with coupling constant g', and Ws^^ are the SU(2) weak isospin gauge fields, with 
coupling constant g. All matrices are the identity in generation space, except Y^^a^ the 
identity in color space, except Aa, and the identity in Dirac space, except 7^. 

The gauge fields shown are the canonical ones. For convenience we shall absorb the 
couplings in the fields, and write the Dirac operator in the form 



in 



SM,q 



ilp^-(f! -^mI 

V i^,-^ -fM] / 



(3.4) 



In this expression v is the vacuum expectation value of the Higgs field after spontaneous 
symmetry breaking, and M^^^ the complex quark mass matrices. 



(3.5) 



Also 



1 



2 cos Ow^^^ 



(3.6) 



(3.7) 



where Z^{x) is the canonical field of the boson and 9w the weak angle. In addition, we 
have introduced the covariant derivatives 



with 



A 



2. 



1/6 



n9 B/j., Ad^n = --gBn 



(3.8) 
(3.9) 



For leptons 



(3.10) 



The lepton field l{x) belongs to the space Wlr Tiue ® "^gcn ® "^Dirac- For convenience 
it includes a spurious right-handed neutrino to achieve greater similarity with the quark 
case. In matrix form the fields are organized as follows 



l{x) 



, l{x) = (^i^L,eL,i'R,eRj , 



(3.11) 
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and the Dirac operator takes the form (we assume massless neutrinos throughout) 



«DsM,l 



-^^Fe 



2^ V3 ^ 






(3.12) 



i ^ 

The right-handed neutrino is completely decoupled. As for the quarks, we find it convenient 
to rewrite the same matrix as 

/O - f -\^^ \ 

-|Me ^JPe + f 


VO -^Ml) 

where we have introduced the new covariant derivative 



^D 



SM,1 



(3.13) 



du + iAe 



For subsequent use, we introduce the derivatives 



(3.14) 



w, 



W-. = Da^^W- - W~Du,f, = [d^, W-] - i(A„,^ - Aa^^)W-, 



f: 



/if 

u,d 



i[D^/,Di:^''] = [d^,A:^'']-[d,,A^/]. 



(3.15) 



Note that the fields are not antisymmetric in fx, v. Also note that the analogous 

construction in the lepton sector gives exactly the same result as for the quark sector, 
namely, = [9^, VF^] ± ig'B^Wj^. In fact, using the relation 



^e.m. 



g'Bn + 2 sin^ 9w (photon field) , 



one finds 



W, 



± 



± 



(3.16) 



(3.17) 



which is covariant under Ue.in.(l)) the remaining gauge freedom in the unitary gauge, apart 
from SUcoior(3). 

3.2 Euclidean space 



We apply to the SM fields the prescriptions given in section 2.2 to go from Minkowskian 
to Euclidean space, and this yields 



D 



SM,q 



D 



SM,1 





V ^, + ^ 

/O ^ + .^ lf+ \ 

|Me W~ 

^0 
Vo % 







(3.18) 







(3.19) 
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with 



Also 



Du,^ = + Au,^,, Dd,^, = + Ad,^,, De,^, = + Ae,f, . (3.20) 



W+ = Du,^W+ - W+Dd,^, = Dd,^W- - W-Du,^, 

fiu 



Fl^l^ = [Dl^\Dl^% (3.21) 



The c-number fields Z^, ^u,/^, Ad^^, A^^^ and F^^ are purely imaginary, while (W^)* = 
—W~ and (VF^)* = —W~^,. The matrix field is antihermitian. The field (/> is a real 
c-number. 



In what follows we shall work in Euclidean space, until section IC, where we return to 
Minkowskian space to display the results. 

3.3 ud-parity 

In the SM only two of the four fields Z^, A^^fj,, and Ag,^ are independent (and F^^, is 
proportional to FJ^^^), however, it proves useful to carry out the calculation for the "extended 
model" where Z^, Au^fj., Ad^^, and A^^^ are independent fields, and also Mu and Md are 
regarded as (x-independent) variables, as this procedure provides simpler expressions. The 
quark sector with generic A^ and Ad enjoys a Uu(l) x 11^(1) symmetry which in the SM 
reduces to Ue.m. (!)• Moreover, for the extended model a symmetry becomes apparent 
under the exchange of labels u and d in the quark sector, namely, 

Mu^Md, W^^W^, Au^Ad, Z^-Z, G^G, (P ^ cj). (3.22) 



This corresponds to a similarity transformation of DsM,q as given in ( 3.18| ) and so it leaves 
the (quark sector) effective action unchanged when expressed in terms of generic Mu, Md, 
Au, Ad and Z (as well a and (p). This symmetry, which we shall call ud-parity, is not 
supported by the SM but it will be present in our calculation and this will become useful 
later. 



4. CP violation 



As noted, the full CP transformation in ( 2.16| ) is a symmetry of the effective action func- 



tional. It is instructive to see this in detail. First, note that, in Euclidean space and in four 
dimensions, the definition of 75 does not contain an imaginary unit i, and so no complex 
numbers are generated in the functional T after taking Dirac traces (nor there are any 
other i's in the Dirac operator or the definition of T, cf. ( p. 14 )). As a consequence, when 



the background fields are replaced by their complex conjugated, T also becomes complex 
conjugated. That is, the (real) normal parity component, is unchanged whereas the 
(purely imaginary) abnormal parity component, T~ changes to —T~. On the other hand, 
the transformation involving {x^,x) — > (x",— x) leaves invariant since it does not con- 
tain e^j^yap-, but changes the sign of T~ . In this way the complete effective action is left 
invariant under full CP transformations. 
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Consider now the (Euclidean) physical CP transformation in the SM. The transforma- 
tion acts on the dynamical fields, namely, Ga, 4>i W^, Z, and Au^d,e, but the constants 
Mu,d,e are unchanged and in general this will no leave the effective action invariant. This 
allows to classify the contributions to the effective action in two types according to whether 
they are even or odd under physical CP: 

r = r+ + r_. (4.i) 



(This classification is not to be confused with the separation ( p.sp into terms even and odd 
under pseudo-parity.) 

In view of the fact that the full CP transformation leaves F invariant, it follows that the 
physical CP transformation has the same effect on the effective action as the transformation 
{v is real): 

M„,,,e-M;rf,,. (4.2) 

This implies the well known result that if the Yukawa couplings were real, or equivalent to 
real, there would be no CP violation. This is automatically the case in the lepton sector 
but not in the quark sector for three or more generations Indeed, for arbitrary complex 
matrices and one can write 

= A~lmuAu,R , Md = AalmdAd,R , (4.3) 

where the matrices m^, nid are diagonal and non negative and Au,l, Au,r, Ad,L, Ad,R are 
unitary, all of them in generation space. Using the freedom to rotate the quark fields in 
generation space allows to bring the Dirac operator to the form 



D 



SM,q 



Inirndfls W ]Pd-^+<^ 
+ ^Q^^muCn^' 



(4.4) 



where 



C = Au,lA^I (4.5) 



is the Cabibbo-Kobayashi-Maskawa (CKM) matrix, and ^^1,2,3 sue arbitrary unitary ma- 
trices (in generation space). It is manifest, by means of an appropriate choice of ^1^2,3, 
that the physics is invariant under the redefinition C U1CU2 where are arbitrary 
unitary and diagonal matrices (so that they commute with and md). On the other 
hand, choosing ^^1,2,3 as the identity matrix simplifies the Dirac operator and shows that 
C has to be complex to allow violation of physical CP; the CP transformation as given in 
( |4.2| ) becomes equivalent to 

C^C*. (4.6) 

The similar manipulations in the lepton sector do not give rise to a matrix C since the 
would-be Ml, complex mass matrix vanishes. So the lepton sector does not contribute to 
CP violation in the effective action and this sector will be disregarded in what follows. 
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Alternatively, the quark rotations in generation space can be chosen so that the Dirac 
operator becomes 



D 



SM,q 



(4.7) 



This shows that only contributions involving can appear in the CP violating sector, 

r_. 

Because all quantities in the Dirac operator are c-numbers in generation space except 
the Yukawa couplings, and the effective action adds one-loop graphs with quarks running 
on them, this functional can be arranged in the form 



r = ^ti fx{Mu,Md) j d^xtiOxix), (4. 



where f\{Mu,Md) are operators in generation space constructed with the complex mass 
matrices and 0\ are local operators constructed with the Higgs and the various gauge fields. 
In the first case the trace refers to generation space, Wgenj and it refers to all the other 
spaces in the second case. As noted before all these operators do not involve any complex 
number in their construction in terms of the fields and mass matrices. As a consequence, 
the CP transformation ( [4.2| ) implies 

tr/A(M„,Mrf) ^ tr/A(M:,M^*) = {iv h{Mu,Md))\ (4.9) 

therefore the CP violating component of the effective action can be expressed as 

r_ = J^iIm(tr/A(M„,Mrf)) j /xtrOA(x). (4.10) 
A •' 

It also follows that the local operators 0{x) contributing to P^t (normal parity, CP violat- 
ing) are antihermitian, since P"*" is real in Euclidean space. On the other hand the local 
operators in PI (abnormal parity, CP violating) must be hermitian. In this reasoning we 
use the fact that the two transformations involved, namely, CP (which defines the separa- 



tion in (4J.)) and pseudo-parity (which defines the separation in ( |2.8| )) commute. This is 
correct since pseudo-parity just exchanges the labels L and R while CP does not mix those 
labels. 

An important remark is that in the Standard Model all bosons can be assigned natural 
parity, (—1)'^. Therefore the abnormal parity sector is just the P odd sector, while the 
normal parity sector of the CP odd component is C odd and P even. This is unlike the 
chirally broken phase of QCD where parity is preserved even by the abnormal parity sector 
due to the presence of abnormal parity hadrons. 
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5. Momentum integrals 



As we shall see, for the operators ti fx{Mu, M^) depending on the complex mass matrices, 
cf. (^^), the derivative expansion produces the set of integrals 



rk 

-^ri,ti,...,rn,tn 



(2^) 



1 



1 



(2^ 



1 „ 1 „ 1 1 „ 1 



C7- 



c- 



(j)2 + m^)''" (j)2 _|_ m^)*" 



(5.1) 



where the exponents k and rj, tj are non negative integers. The integral with 2n indices will 
appear in contributions with n and n . (Charge conservation requires as many 
as W~ in any contribution to the effective action and are the only fields connecting 
the spaces u and d.) 

Due to the cyclic property of the trace 



and also taking, the hermitian adjoint. 



/jfc )* = j'^ 



(5.2) 



(5.3) 



Since we have seen above that the CP violating component is tied to the imaginary 
part of this integral, we introduce the definition 



fk 

^r\,ti,...,rn,t„ 



^-'■™-^ri,ti,...,r„,t„) 



which enjoys the properties 



'ri,ti,...,r„,t„ 



r2,t2,---,1'n,tn,,ri,t\ 



;^n— 1 — 1 V" 5^1 5^1 



(5.4) 



(5.5) 



From these relations it is immediate that I vanishes if n = or n = 1. Therefore, 
at least 2 C and 2 are needed to have a contribution to r_, or equivalently, 2 
and 2 . This is a well known fact in the literature Q. (This implies that the operator 



tr [F^yF^'^) [20 1 mentioned in the Introduction cannot be derived from simple integration 
of the quarks as the would-be term with four vanishes. Such term can be produced if 
internal gauge field lines are allowed.) 

Of particular interest will be the first non trivial case, n = 2. For it one finds 



ri,ti,r2,t2 



r2,ti,ri,t2 



ri,t2,r2,ti ■ 



(5.6) 



It will also be useful to note the transformation of the momentum integral under 



tid-parity (see section |3.3|) , namely. 



Tk 

^ri,ti,...,rri,t„ 



rk 

'^ti,...,r„,tn,ri ! 



(5.7) 
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and in particular 



^ri,ti,r2,t2 ^ ^h,ri,t2,r2 ■ (^•^) 



Using the second form in (5^), to compute the integrals of the type /^^ ^^ ^^ t2 three 



generations, one can use the identity |^] 

Im {CijC~i}CkiCf^) = Jeikeji, k,l = 1,2,3, (5.9) 

(with no implicit summation over repeated indices here) where €ij = J2k=i ^ijk, and J is 
the Jarlskog invariant [P, |l^] 

J = cos 012 cos^ 013 COS 023 sin 012 sin 0i3 sin 023 sin 5 = 3.0(2) X 10"^ (5.10) 

The resulting momentum integrals no longer involve matrices in the integrand. 
The only integral required in this paper is /f i 2 2- It can be cast in the form 

ili^2,2 = iJGFKcv ■ (5.11) 

KCP is a dimensionless coefficient and Gj? is the Fermi constant, which can be related to 
the the vacuum expectation value of the Higgs field, Gp = l/(\/2f^). For v and for the 
quark masses we take v = 246 GeV, niu = 2.55 MeV, = 5.04 MeV, rus = 104 MeV, 
rric = 1.27 GeV, rrife = 4.2 GeV, mt = 171.2 GeV 0. This gives 

Kcp = 3.1xl0^ Kcp = 3.6 X 10"^ GeV"^ (5.12) 

The integral 1 2 2 homogeneous function of the quark masses of degree —2 so 

KCP can be expressed as a function of the Yukawa couplings, Uq = ^/2mq/v. The loop 
function kcp is of great interest by itself so we give full details of its form and calculation 



in section |11|. We only remark here that it is numerically many orders of magnitude larger 
than the "minimal" term {yl-yl){yl-yt){yt -yl){yl-y1){yl-yl){vl-yl) = 6x 10"^^. 
The reason, of course, is that other non minimal factors appear in the full expression. 

6. Derivative expansion in the SM 

Smit [|^] has proposed to use the derivative expansion as a suitable approach in the present 

context of CP violation in the SM. 

In the SM each , Z , D^^d,!/ or Ga counts as first order. As we have seen, at least four 
are needed in the CP violating sector, therefore r_ vanishes at zeroth or second order 

in the derivative expansion. It was shown in [11| that the fourth order is also vanishing in 

the abnormal parity sector. We verify below that the fourth order vanishes actually in both 



sectors. Therefore, as suggested in [11|, the first non trivial contribution should start at 



six derivatives. Hernandez et al. |12] have addressed such a computation in the abnormal 



parity sector and find a non vanishing result. We compute below all contributions to sixth 



order, including Higgs and normal parity terms. Our result do not sustain those in |12|. 
We find non vanishing contributions in the normal parity sector but none in the abnormal 
parity one. 
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As we have noted above we do not need to consider leptons since they do not give a 
contribution to r_ assuming massless neutrinos. On the other hand we can also neglect 
gluons. At fourth order the four already saturate the required number of derivatives 
and no gluons are allowed. At sixth order there is room for up to two gluon fields. However, 
by gauge invariance one gluon is not admissible and two gluons must be combined to form 
a gluon field strength. Such term vanishes due to the trace on color. In what follows gluons 
are not included and color just gives a global factor Nc = 3. 

Another question is whether the derivative expansion at low orders produces a reliable 
approximation to the physical amplitudes. Using simple estimates, it has been argued in 
1 12] that the range of validity could reach the scale of the charm quark mass or even larger. 
Besides, it is clearly of interest to have correctly accounted for the lowest order operators 
of the effective action, as guidance on the available CP violating mechanisms. 



7. Chiral invariant approach to the effective action 



The approach of section 2.2 can be directly applied to the SM. 

For the SM in the quark sector, the operator K of ( p. 19 ) takes the form (the gluons 
are no longer present in the covariant derivatives) 



Here we have introduced the shorthand notation 

^^{x) = cl>-^[d^,4>]. (7.2) 

The operator K acts in the space TCud <^ "^gen ® '^coior ® '^Dirac (that is, the same space as 
DsM,q except the factor Hlr)- Therefore, for the SM the equations (|2.20| ) become 

r+ = -iiV,Re(logK), 

r- = -^N,ilm{logK)+T^wzw- (7.3) 

Here we have introduced the symbol ( ) which will be used in what follows. It denotes 
a trace operation including a 75 in the abnormal parity sector, and just the trace, without 
75, in the normal parity sector. The precise trace operation implied by ( ) will often be 
obvious from the context. The inclusion of 75 does not spoil the cyclic property for ( ) 
since all operators involved will have an even number of Dirac gamma matrices. In this 
way we can treat simultaneously the normal and abnormal parity components. 
In ( |7.3|) the trace implied by ( ) acts on x-space and on Tiud ® Tlgen ® "^Dirac- 
It has been shown in [^] that the term T^y/zw does not have a contribution to the 
CP violating component of the effective action, so we disregard this term in what follows. 
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For convenience, let us separate K into its diagonal and off diagonal parts (in ud space) 



K = Kd + Ka, 

( {<pyv^)MuMl - {Ip^ + ^{p^ + jp) 



(0V^2)^,m]-(^,-^(^,+ 



The advantage of this separation is that only contains charged currents, which can 
break CP, and also that this term is of first order in derivatives^ (recall that the gauge fields, 
and in particular W^, count as first order) so the n-th order in the derivative expansion 
of r can contain at most n factors K^. 

Substituting this form in (logET), which appears in (^^), yields 

VJ-R-.J-R-.\-i/.i ^-1^-1-^1 



(logK) = (log A-o) - - [^K,-K,J - - (--K,-,U-K,-K, 

00 

= ^(logK)2„. (7.5) 



n=0 



The subindex 2n in (log-fr)2„ indicates that the term contains exactly 2n fields. Work- 
ing out the trace in ud space (for n > 0) one obtains 



(7.6) 



Here the trace implied by ( ) refers to TL^en ® Woirac, and x-space. This expression is 
manifestly ud-parity invariant (see section ^^), as can be verified by using the trace cyclic 
property. 

Because at least four are required to have a CP violating term, the relevant 
contributions start at (log 7^)4. To sixth order in the derivative expansion only (logi^)4 
and (log 1^)0 have to be retained. (log-fr)4 contains terms with at least four derivatives 
(namely, those coming from M^^), while (log-ftr)g starts at six derivatives: 

(log 7^)4 = (logi^)4+o + (log/04+2 + 0{D^), 

(logi^)6 = {\ogK)^^, + 0{D^). (7.7) 

Here (logi^)2„+2m indicates 2n derivatives from and 2m more derivatives not from 
W^^ i.e., coming from 9^, or A^''^ . 



^Ka is of order one and not of order two because the counting refers to derivatives of the external fields 
and in (^^) Ip^ ^ can still act on the quarks. As noted previously in section 2.5, such derivatives on the 
running fermion are of order zero in the derivative counting. 
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8. Vanishing of terms of the type 2n + 



In Ref. |11|, and also confirmed in |12|, it was shown that there is no CP violating contri- 
bution to four derivatives in the abnormal parity sector. Let us show that this is true in 
both sectors and that, moreover, there is no sixth order contribution either coming from 
terms with six W^. 

These statements are remarkably easy to establish using the method of symbols de- 
scribed in section 2.5. This method amounts to make the replacement d^^u d + |^ in 
Kd and Ka and integrate over p^. (Recall that is purely imaginary but p"^ = — p^, as 



explained in section |2.5| .) 

Concretely, the CP violating terms with precisely four derivatives must come from 
(logi^)4 taking no other derivatives than the four (i.e., must be of the type 4 -|- 0). 
Therefore, to four derivatives, we can set Ip^ d^i^-- *-* f/p-^ in the operator K: 

1 f d'^xd'^p 



(log i^) 4+0 



2 J {2ttY 



(8.1) 



Here it is already obvious that, upon momentum integration, the integral 1 1 1 (in- 
troduced in (5T)) will be generated. Because this integral is real, due to eq. (|5.3D, it 
follows (cf. (|4.1C| )) that no CP violating term is produced to fourth order in the derivative 
expansion, neither in the normal nor the abnormal parity sectors. 

To see this in more detail, we first take an angular average in (^.1[), using 



ViiVvVaVp {^iiv^afi + ^iia^vji + ^^iP^au)?^/ {d{d + 2)). (8.2) 

Since no derivatives with respect to x are present in the expression, we can simply rescale 
P^l {4'ix)/v)p^, and the momentum integrals in (5.1) apply. Specifically, 

(logi^)4+o = I d^x { (r+7M W-l^.f + [W^1^. W~l.f 

+ W'lu w^iu r"7M)- (8-3) 

The result to four derivatives is proportional to ll as advertised. 

If we consider now the case of sixth order 6-1-0, i.e., when the six derivatives are 
saturated by six , it is quite clear, by using the same reasoning, that the result will be 
proportional to if ^ m which is also real, and therefore, also no CP violation is produced 
in either sector from such contributions. Of course, the analogous result holds for all orders 
of the type 2n -|- too. That is, 

(r_)2„+o=0. (8.4) 

We note that this vanishing is rather trivial in the abnormal parity sector (for space- 
time dimension d> 2): with the Levi-Civita tensor and only two four- vectors it is not 
possible to construct a non vanishing scalar. 

The vanishing for the normal parity part is also easily understood. Due to charge 
conservation, the possible operators constructed using only are of the type 
{{W+W+){W-W-)Y{W+W-)'^ and are CP even. 
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9. CP violation in the absence of Higgs field derivatives 



We have just seen that no CP violation occurs to four derivatives and also to six derivatives 
if these are saturated by T^^'s. Therefore, any CP violation through sixth order in the 
derivative expansion must come from terms with four and two other derivatives not 
of the type, that is, terms of the type 4 + 2: 

r_ = (r_)4+2 + o(i?'). (9.1) 

In the next section we shall consider the general case. Presently, we study the simplest 
situation where no derivatives of the Higgs field are considered. In this case, the Higgs 
field (f){x) itself can be set to its vacuum expectation value v, since it can be restored in 
the formulas at the end by a rescaling of the quark masses. Under these assumptions, the 
trace (logi^)4 of (7.6) reduces to the simpler form 

{logK), = (M„Mt - (^„ + ^ p^y' W^JPd {m,mI - (% - ^ p,y' W'Pu 

(9.2) 

Here we apply once again the method of symbols making the replacement — > 
+ Dfj^, integrating over and then expanding in powers of the derivatives: 

(9.3) 

Where 

= {M.,Mi -[i>+ip^ + m + jpu)y\ = {MaMi -{i> + ip^- m + p^t'- 

(9.4) 

Expansion of the first denominator gives 

(9.5) 

and similarly for N^- And we have defined 

Nu = {MuMI+pY\ Nd = {MdM^+pY^ . (9.6) 

The two objects A^^ and are 2;-independent, they do not commute with each other but 
commute with all other quantities in (logK)^. In addition, they appear in the momentum 
integrals introduced in section ^ Indeed, the first eq. in ( ^.1| ) can be rewritten as 

At this point considerable simplification can be achieved by making the following 
observation. To produce (logEr)4_,_2 from (|9.3| ) we need to pick up exactly two derivatives 
(apart from the four explicit PF^). On the other hand, CP violating contributions come 
only from the imaginary part of /^^ t2' ^^'^ ^^^^ requires ri ^ r2 and ti ^ t2, cf. (^). 
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Now, it is quite clear that, in order to obtain such a situation, it is necessary to pick up 
exactly one of the derivatives from one of the A'^ and the other derivative from one of the 
Nd- Any other possibility ends up with either the two A^^ or the two raised to the same 
power, i.e., ri = r2 = 1 or ti = t2 = 1- Therefore the CP violating contributions in the 
present case come from /f i 2,2- 

To work this out let us simplify the expressions by introducing the following quantities, 
which appear naturally in (^]^) and ( |9.5D 



w^=\^^, 5u = 2pDu+^i>, 5d = 2pDd- fTt> . (9.8) 
Then, applying the previous observation yields 

(logi^)4+2 = "2/ -^^{^I^uw+Nlddw-N^w+Ndw" +Nl5uW+NdW-NuW+Nl5dW~ 
+ CP invariant terms 

X (^5uW^ 5dW^ — 5uW^ 5dW^ — w^ddW~5uW^w^ 

+ w+w-6uw+6dw-'^ + CP i.t. (9.9) 

In the second equality we have rearrang ed the factors N^, N^, Nd, N^, using that the CP 
violating part of the momentum integral, I^^ ) antisymmetric under exchange of 

ri,r2 or ti,t2- 

The integrand in (|9.9D contains derivatives (inside 6u,d) which are not derivating any- 



thing yet. As explained in section 2.5, in general one proceeds by moving the derivatives 
to the right, and at the end the momentum integral kills these "free" derivatives. In the 
present case this turns out not to be necessary. Instead, we can introduce the combinations 

(dw)^ = 6uW^ — w^Sd, {6w)~ = 6dW~ — w~6u, (9.10) 

in such a way that 

(logK)4+2 = "2/ -^triNuNdN^Ni){i6w)+w-i6w)+w- -w+i5w)-w+iSw)-) 
+ CP i.t. (9.11) 

and the integrand no longer contains any "free" derivative. 

The CP violating part of this result is proportional to if 122 develops a factor 
{v/(f){x))'^ upon restoration of the Higgs. 

There is a simple observation that can already be made at the present stage. Namely, 
there is no CP violating contribution to the abnormal parity sector from terms of the type 
4 + 2 without Higgs field derivatives. As discussed in section 0, in the abnormal parity 
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sector the operator multiplying the momentum integral must be hermitian to have a CP 
violating contribution, however, 



O = tr [75 {{5w)+w-{6w)+w- - w+ {5w)- w+ {5w)-)] (9.12) 

is purely imaginary. To verify this, we take the complex conjugate of everything inside O. 
In Euclidean space four-dimensional space 7^ and 75 are related to their complex conjugates 
by a common similarity transformation, 

i; = C-'j^C,, 7* = C-S5C,. (9.13) 

Also, one verifies that 

{w^y = C^^w^Cc, i{5w)^)* = -C-\5w)^Cc. (9.14) 



Therefore 



as advertised. 



O* = -O, (9.15) 



In the absence of Fji^'^ and of complex quark mass matrices, complex conjugation 
becomes equivalent to ud-parity.'' So O is imaginary (in the normal and in the abnormal 
parity sectors) because it is odd under ud-parity. In turn this was obvious without further 
calculation once the momentum integral 1 2 2 obtained in (|9.9|). This is because the 
latter is odd under ud-parity and the full effective action is even (cf. section |3.3D . 

Remarkably, the operator tr [j^{5w)^w^ {Sw)^w^] (and hence its complex conjugate) 
vanishes by itself after taking an angular average and the Dirac trace. We have not found 
a simple explanation for this. 

The operator O in the normal parity sector (i.e., as in ( |9.12 ) without 75) is also odd 



under nd-parity and so also purely imaginary. Therefore the operation of taking the real 
part indicated in ( |7.3| ) is redundant. The normal parity contribution is not vanishing. The 
result so obtained is part of the general result which we present in the next section. 

10. CP violating terms to six derivatives 

In this section we present the full result for the CP violating terms of the effective action 
to six derivatives. This includes all relevant fields in the SM, and derivatives of the Higgs 
field. 

We have used the method of symbols and repeated the calculation using the method of 
covariant symbols as a check, to obtain precisely the same result from both calculations. In 
the latter case we use the covariant derivatives Du,^ and -D^,^ to carry out the construction 
indicated in ( 2.26| ). This full result is also consistent with the independent computation 



made in the previous section. 

From the calculation we obtain the remarkable result that {logK)^^2 vanishes iden- 
tically in the abnormal parity sector for all terms that could have a contribution to CP 



''Under complex conjugation W,^ —W^, Zf_i ^Zf^ (in Euclidean space) while under ud-parity 
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violation. In fact, T~ vanishes for all the terms that we have computed, whether CP vi- 
olating or not. (Of course we have not studied most CP invariant terms with six or less 
derivatives, as they are not required for our purposes.) We have not found a compelling 
reason for this, so most likely, the vanishing found is just a low order accidental symmetry. 
The existence of CP violating terms in the abnormal parity sector of the SM with eight 
derivatives or more is not excluded. These would be the leading P violating contributions, 
relevant to the electric dipole moment problem. 

Another unexpected result is that {\ogK)^_^_2 is purely real in the normal parity sector, 
for terms that contribute to CP violation. (And so, taking the real part indicated in ( [7.3| ) 
becomes redundant.) In the calculation this follows from the fact that only the momentum 
integral if i 2 2 appears. This integral is odd under ud-parity and hence the accompanying 
operator must be odd too and hence imaginary. Once again it not obvious to us whether 
this feature will be maintained at higher orders in the derivative expansion, coming from 
an exact selection rule in the SM, or is just an accidental symmetry. 

The result, in Minkowski space and in the unitary gauge, reads 

TsM = -^iJ Gf Kcp j (fx (00 + 01 + 02) +C'(-D^) + CP invariant terms. (10.1) 

Here A'^c = 3 is the number of colors, J the Jarlskog invariant, Gp the Fermi constant and 
= 3.1 X 10^ is the dimensionless parameter of section ^. The operators Oj, i = 0, 1, 2, 
have dimension six and are all purely imaginary: 



2 2 

On = -W^W~W^W~ W^W W^W 

4 

J — W^W^ W^W~ - 2W~^W~ W^W~ 

2 1 
K 1 

+^W+W-JV;^W-^ - -W;W-^W;W-, - c.c. 
Oi = ^(v^^ + iZ^) 

■^^ \" L/ " fj, " u " aa " u " u " fi " aa 

+ W^W~'W^W~ - W^W'W^W" 

I '"i/"a"^"au " h> " a " u " ^a 

- W+W-W+W~^ + W+W-W+W-J) - c.c. 

-^{ip^ + iZ^){ip^ + iZ^) (10.2) 
X (W+W-W+W- - 2W+W-W+W- + 2W+W-W+W-') - c.c. 

In these expressions c.c stands for complex conjugate. Even if these expressions refer to 
Minkowski space, the Lorentz indices are all written as subindices for clarity as no ambiguity 
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may arise. With the conventions given above, the pass from Euchdean to Minkowskian 
metric amounts to the replacement — > iZ^ with no other change. 

As noted in section 2.4, the effective action in Minkowski space is purely real, at every 
order in the derivative expansion and this property is found here. 

The fields were defined in section 3^ and they are expressed in terms of the 
Ue.m.(l) covariant derivative in ( p. 17 ). On the other hand ^p^ was defined in ( [7.2[ ) as the 
logarithmic derivative of the Higgs field. 

As it is readily verified, all the operators Oo,i,2 are indeed odd under the CP transfor- 
mation 



W^{x) ^ -vr/ W^{x), W^,{x) ^ -^/^/ W^pix) . (10.3) 

Terms including the field strengths Fji^"^ are absent. This can be understood from the 
fact that the available operators, iFji^'^W^W'W^W' , are CP even. 

On the other hand, the result in ( |10.2| ) presents some regularities which for us remain 
purely "empirical". In Oq terms coupling W^^ to do not appear. It is always possible 
to change variables from cp^ and to ipf^ztiZ^, however, it is not obvious why, in Oi, the 
combinations (/?^ ^iZ^ couple only to and not to W^^. Also it is not clear why in O2, 
the combinations 93^ it iZ^ couple only with themselves and not with (pfj_ =p iZ^. This latter 
observation suggests the speculation that the effective action (or perhaps T~^) to all orders 
could be of the form F[ip^ + iZ^] — F[ip^ — iZ^], where the functional F would depend 
analytically (holomorphically) on its argument. 

Also, the relatively simple dependence of the result on the combinations 99^ it iZ^ 
suggests the possibility of reconstructing the full result with Higgs derivatives from that 
without (p^ (by some kind of gauging). In this case the calculation in section ^ could 
perhaps be adapted to include ip^. We have not tried this in this work. 

The fields (p^ — iZ^ and iW^ follow from projection of V^^> onto $ and <I>, respectively, 
where represents the full SUl(2) x Uy(l) covariant derivative, and W^^ can also be 
written using ^ and <I> and their covariant derivatives, but the result is not particularly 
illuminating. 



11. The coefficient kcp 

In this section we study in some depth the function kcp- Using the second form in ( [5.1| ) 
as well as the identity ( ]5.9| ) , the momentum integral -/^i 1 2 2 takes the form 

*3 _ ■ T ^ /" d'^'P . 9^^^ 1 1 1 1 



1,1,2,2 ^ik^jl I (o^M^^ ' {^2 \ r^2 _L ^2 ) (^2 _^ ^2 ^)2 (p2 ^ j„2 ^)2- 

(11.1) 

Here mu,i denotes the mass of the quark of type u of the i-th generation, and similarly for 
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The sum over flavors is easily carried out using the identity 



y^^^ 1 1 _ {ml - ml){ml - mf){m^ - ml) 

* (p2 + ml j) (p2 + ml j^y (p2 + m2)2(p2 _|_ ^2y(^p2 _|_ ^2-j2 ' 

and similarly for d-type quarks. The integral can then be written as 

-^1,1,2,2 ~ ^mlm- (11-3) 

where 

Am = {ml - ml){ml - mf){mf - ml){ml - m^)(m^ - ml){ml - ml) , (11.4) 

The factor J Am is just the Jarlskog determinant A of ( |1.1D . This is the "minimal" fac- 
tor that must always be present in the CP odd effective action as dictated by perturbation 
theory. This can be seen for instance considering the same integral at finite temperature 
which amounts to replace the energy integral by a fermionic Matsubara sum. In the limit 
of sufficiently high temperature all quark masses can be treated as perturbations. So the 
known perturbative result applies and it starts at order twelve with JA^- As a conse- 
quence the integral if i 2 2 consistent with the well known fact that no CP breaking can 
take place if two up-type or two down-type quark have the same finite mass. (The case of 
degeneracy with vanishing mass requires a separate study since Im could present infrared 
divergences.) 

In Im the product is over the six quark flavors, regardless of its u or d type. So the 
coefficient kqp introduced in section |5| has the same symmetry as A^ under exchange of 
quark flavors, the factor Im being completely symmetric, and actually positive definite. 

The momentum integrals of the type Im, typical of zero momentum insertions in a 
Feynman graph, can be computed using the relation 



where 



J (2^^^'^' n (p2 + „2)r, = (4vr)'^/2r((i/2) A; + d/2 = 1, 2, 3, . . . 

J — 1 J 

(11.6) 

I-,...,r„{mu . . . , m.) = ^ log(z) n , (11-7) 

j=i ^ j' 

and the integration is along a positive closed simple contour enclosing the poles at m| but 
excluding z = 0, and the cut is on the real negative axis. The identity (jlL|) assumes 



positive m| and holds whenever the left-hand side is ultraviolet and infrared finite. If it is 
not, the right-hand side gives the finite part.^ 



*Note that when a is a non negative integer, the contour integral scales as 1 + a — "Y^- rj but in the 
presence of infrared divergencies an anomalous scale term develops from the logarithm . 



-25- 



Using ( |11.6D , and comparing with the definition of the loop function kcp in ( 5.11] ) 
yields 



KCP = 2^/2 Ay ly, (11. 



where 



A, = {yl - vlM - vhivl - ylM - vDivl - yl){yl - vl) = iV2/vr a^, 
= j^:^^2,2,2,2,2,2{yu,yc,yt,yd,ys,yb) = iv/V2)^'^im, (ii-9) 

and the Yukawa couplings yq = \/^mq/v are used. The coupling controlling the strength 
of the CP violating operators has dimensions of one over mass squared. The use of the 
Yukawa coupling amounts to using v as convenient mass scale to measure this coupling. 



The contour integrals ( 11.7 ) are readily computed by residues.^ This produces the 
explicit expression 



yfjl + Uogyf) _ yflogyf / § _ o, 6 V-/ 



yl 



+Qyf E' y'yl - ^y" E' yUy" + ^ E' yUy'yr 

j<k j<k<l j<k<l<r 



(11.10) 



In this expression the indices i,j, k, I, r run over the six quark flavors. The prime in a sum 
or product indicates to omit the term i in that sum or product. Although the explicit 
expression of the integral ly looks divergent as two masses become degenerated this is 
not so, as is obvious from the integral itself: in the coincidence limit the integral ( 11. 5| ) 



is perfectly regular. Also, ly yields an homogeneous function of the y of degree —14. 
The inhomogeneous scale variation from the logarithms cancels after adding the six terms. 
(This implies that one can use one of the Yukawa couplings, yq, as overall scale, yi — > yi/yq 
to remove one of the six logarithms in the expression, at the price of a less symmetric 
formula.) The formula can be written in many different ways but none is expected to give 
a simple expression. The simplest and most transparent form is perhaps its very definition 
as a momentum integral, ( |11.5| ). 

Each of the factors Ay, ly in kcp gives quite disparate numbers, namely. Ay = 6.0 x 
10^^^ and ly = 1.8 x 10^^. This is because they are homogeneous functions of very large 
degree (high mass dimension) and so they change wildly under even moderate changes in the 
scale. For instance, in units of the bottom quark mass A^, = 1.5 x 10^ and = 4.1 x 10~^. 
Their product has degree —2 and so the number is much less dependent on the mass scale. 

Therefore the huge cancellation between factors is partially trivial, and as noted in |jll| 
the factor Ay should not be used as a rough estimate of the CP violating component. What 
is not trivial is the detailed role played by the light quarks. The value of kcp is enhanced 



(m^)°log(m?; 
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by the small mass of the quarks m, d and s. This can be seen by changing artificially their 
mass, niu XrUu, rud — > Am^, rUg Xnig- As A moves in the range from 1 to 25 one finds 
that ly decreases monotonously by a factor 5000. This rather large factor corresponds to 
an effective dimension which runs from about —2 at A = 1 to about —2.5 at the highest 
value of A. The value —2 corresponds to the scale dimension at the infrared divergent point 
A = 0. On the other hand the factor Ay increases by a factor about 400. Beyond this range 
Ay reaches a maximum and starts to fall heading for the zero at A = ray^jms ~ 40, and then 
starts to grow again with opposite sign. The net effect is to find a quenching of kcp as one 
departs from the SU(3) chiral point. Nevertheless, it should be noted that kcp is a rather 
complicated function of the quark masses, without well defined sign, so any interpretation 
should be taken with some caution as the result might look different depending on how the 
masses are moved in detail. 

The coefficient kcp is infrared finite as any two quark masses go to zero while the 
other quarks remain massive. This allows us to consider the limiting case m„ = = 0. 
(The coefficient vanishes trivially if the two massless quarks are n-like or d-like, due to the 
antisymmetry of the Aj, factor.) This gives for the ly factor 

..2 



(47r) 



ttV ^ yfUiiyf-yi)^ IlM-yi)^^ , f^k 



Here the indices k run over the four remaining flavors s, c, b, t. 

The approximation m^, does not change much the value of the coefficient kcp, 

overestimating it by a 3% as compared to the exact value. Since the mass of the quark 
s is also rather small one can consider the further limit — ^ 0. The limit of vanishing 
mu,md,ms is affected by infrared divergencies and this manifests in the fact that while 
Kcp remains finite, the value depends on how this limit is taken. Explicitly, the dominant 
term as mu,md,ms become small is 



Ayly 



1 



(47r)2 



-yjvM + y's) + ^ylyhM + vt) - vM + y's) - ylivl - yl? 
{yl-yDiyl-yiy^yl-ylf 



-'^y>z{yllyl)j^-Mr^. (ii-i2) 



+2yflog(y,V2/n 



(2/1 - yifkal - ylf 

, '^yWd - yWs - ylvs 
(yl - ylfiyl - yl? 



(yy2 y2 ^ + ^ (2/Ught quark) • 



The terms retained are those which are homogenous of degree zero in yu,yd,ys, while the 
reminder has degree 2. This function is antisymmetric under the exchange of d and s. A 
quite remarkable fact is that, to this order, the dependence on yc and yt factorizes and 

moreover 1/5 does not appear in the expression. 

The above function is finite but far from continuous in the massless limit. The phys- 
ical situation is m„,md <C so a sensible limit to consider corresponds to taking first 
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n^'ui'iTT'd ~^ and later rus 0. In this approximation 



Equivalently, in this SU(3) chiral limit 



Gf /^CP = - A ) ■ (11-14) 



The full result corresponds rather to 0.92/((47r)^m^) so this approximation overestimates 
it by an 8%. Despite the simplicity of this approximation the massless SU(3) chiral limit 
is highly non trivial as can be seen from ( p,1.12| ), and in fact even the sign is not obvious 
from the expression. 

12. Conclusions 

We have studied the CP breaking terms in the effective action of the Standard Model 
obtained by integration of quarks and leptons, including operators up to dimension six. 



The result of the calculation is summarized in eqs. ( |10.1| ) and ( |10.2D . 

One of the main results is that such CP breaking operators appear with a sizable 
coupling, of the order of 5 x 10^ times the Jarlskog invariant times the Fermi constant. 
This is much larger than predicted from Jarlskog determinant considerations based on 
perturbation theory and is fully in agreement with the expectations first put forward by 
Smit in |11|. 

Remarkably the non vanishing CP violating contributions come from the normal parity 
sector. This is is somewhat unexpected as it is usually taken for granted that the presence 
of the Levi-Civita pseudo-tensor is needed to have CP breaking. Also noteworthy is our 
finding that, to the order studied, all abnormal parity terms vanish. This is in conflict with 
the result presented in |]l2| where a non vanishing abnormal parity contribution is derived. 
Our result implies that the first CP odd and P odd contribution, relevant to electric dipole 
moments, requires at least the eighth order in the derivative expansion. 

The full result presents interesting regularities (including the just mentioned vanishing 
of terms which are simultaneously CP odd and parity odd) which may follow from the 
structure of the SM or may be accidental symmetries surviving only at lowest orders in the 
derivative expansion. 

The calculation presented here applies to zero temperature. At the order studied in 
this work, the coupling is controlled by the integral if 12 2- The same integral can be 
considered at finite temperature replacing the energy integral in the loop by a fermionic 
Matsubara sum. At high enough temperatures such integral becomes proportional to the 
Jarlskog determinant times the accompanying power just by dimensional counting. 

For temperatures comparable to the top quark mass, relevant for baryogenesis, this mass 
can no longer be treated as a perturbation. Nevertheless, the integral becomes equal to 
the Jarlskog determinant times an order of unity function of mt/T and the required power 
l/T^^, so the numerical estimate is similarly small. At higher order in the derivative 
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expansion other momentum integrals will be generated but the analysis is expected to be 
similar. On the other hand, at finite temperature new operators can be produced even at 
the same order studied here. Due the to breaking of Lorentz invariance down to rotational 
invariance [32| new combinations may arise, e.g., with abnormal parity. 



The effective action obtained here can not be applied directly to CP violating hadronic 
processes since the quarks have been integrated out. In this view it would be of interest to 
repeat the calculation but inserting external hadronic currents in the quark sector. This 
would allow to account for CP violating contributions to meson decay amplitudes. CP 
violation requires to switch between the u and d spaces and in the present calculation this 
can only be achieved by the action of which are the only charged particles left after 
integration of fermions. The situation may change in the presence of hadronic insertions 
carrying electric charge, and presumably lower order operators would be allowed replacing 
some charged gauge bosons by hadronic currents. 

Related to the above, we have seen that the coupling kcp controlling the strength of 
the CP violating operators is a complicated function of the quark masses and is highly 
non continuous in the relevant limit of taking the light quarks to be massless. This is 
due to the presence of infrared divergences in the chiral limit. This suggests that the 
standard chiral perturbation theory corrections induced by the pseudo-scalar Goldstone 
bosons could introduce sizable modifications to the result. Formally, we have included QCD 
in our result, although no gluonic correction have been needed to the order considered. Of 
course, the calculation can be organized under different schemes, for instance counting 
independently derivatives coming from electroweak fields and QCD fields, being the latter 
either fundamental or hadronic. A point to note in this regard is that in the presence of 
hadrons abnormal parity is no longer equivalent to parity violating, as the pseudo-scalars, 
for instance, have abnormal parity. 

Likewise it would be of interest to consider the effective action in the sense of the 
Legendre transform (after integration of fermions), i.e., adding one particle irreducible 
graphs. This could produce lower dimensional CP breaking operators, such as the operator 
tr {F^^F'^'^) mentioned in the Introduction. This can be so by integration of some of the 
external lines in the loop (hence making these lines internal) either in a single fermion loop 
having CP violation or by coupling more than one loop, one of them carrying CP violation. 

As a final comment we note that leptons have not been included as no CP breaking 
takes place in that sector when the neutrino masses vanish. However, our calculation 
applies to Dirac neutrinos with non zero masses, and in this case the leptons produce a 
CP violating contribution completely similar to that obtained for quarks, replacing the 
Jarlskog invariant of the CKM matrix with that of the Maki-Nakagawa-Sakata matrix p6| , 
and using lepton masses to produce Kq^'"'^''. This function will be identical to that discussed 
in section 11 for quarks. Similarly as for the case of quarks, the coefficient not 



at all continuous as some of the leptons become massless. This implies that the limit will 
depend crucially on how the small neutrino masses compare with each other. 
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